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INVESTIGATION  OF  POLAR  MOTION  FROM  DOPPLER  TRACKING  OF  THE 
NAVY 'navigation  SATELLITE  SYSTEM  DURING  THE  MERIT  CAMPAIGN 

WILLIAM  H.  WOODEN,  JOHN  A.  BANGERT.  AND  J.  MILO  ROBINSON 

DEFENSE  MAPPING  AGENCY 
WASHINGTON,  DC  20305,  USA 

APR  2  9  1986  - 

summary 

For  more  than  a  dozen  years,  Doppler  tracking  of 

Navy  Navigation  Satellite  System  (NNSS)  has  demonstrated  the  . ' 

ability  to  give  solutions  for  polar  motion  as  a  byproduct  of 
precise  ephemerides  generation.  This  paper  reviews  the 
polar  motion  estimation  process  at  the  Defense  Mapping 
Agency  from  Doppler  measurements  of  the  NNSS  and  gives 
comparisons  to  polar  motion  series  from  newer  observational 
tecltniques  and  to  the  Bureau  International  de  1  Heure  (6IH) 

Circular-D  series.  The  different  techniques  are  compared 
for  the  period  of  the  MERIT  Main  Campaign.  In  addition  to 
comparing  results  for  the  Nova  and  Oscar  NNSS  satellites, 
comparisons  are  done  utilizing  the  WGS  B4  and  the  NWL  922 
gravi tati onal  fields  in  the  precise  ephemerides  generation 
process . 


1 . 0  INTRQD'JCTTQN 

Since  the  late  nineteenth  century,  it  has  been 

recognized  that  the  i nst antaneous  ^pin  axis  of  the  Earth 
moves  with  respect  to  the  geographic  pole  of  the  Earth's 
crust.  This  "polar  motion"  was  predicted  by  Euler  in  1752, 
but  was  not  conclusively  observed  until  Kuenstner  s  work  in 
The  International  Latitude  Seirvice  vILSi  was 
er.:  c  ab  (  1  in  to  cont  i  fiuous  I  v'  monitor  the  morion  of 

n  sysi  eif'-- •  !•:  detar^m  nat  i  oi  is  at  latitude. 

Un  tt'  the  •le-.’v?  I t  n~-  .o-w  ^  st‘-o;:u»t  r  r  i  fi  s  c’-iun-ai^t  s  during 
til*-"  !.  "ii  ar'iii  use  a:  tna  n-'.ioiia-.  onser  v- tor  i  es .  it 

W.-I  -  df/n::  I.  lie' I  t ''i  r ei" '■■cnjn  i,  .’a  >.  I'u?  i  LE  im.o  Lira  T •  i  ter nat i  onal 

pill  -r  >1'in  I.  tr  r  i  ni.j  ofr-'vioi-.*  (  i.  .-'ii''- >  in  1^32  wi  t.r  ■  ne  charter  of 
der  i,  /iriij  f.io'rir  i  on  +roir'  latituvJe  arid  uroversal  time  data 

(•  1 1  a  L  I.  A  t  r  on  am  :  a  i  i  n  s  t  r  lj  irien  t  s . 

Tne  Hur  ear  I r.  t r-'r  riat  i  onal  de  I  Heure  '  Bii-l.'  was  created 

;.n  1. 1 2  to  nnr*,'  t  ..  me  b'/  putlishm-j  the  time  o emssion  ',in 

r.i  I.  -  I- 'real  i  i'lio,'  o-:-  r  a'ir  o  t  i.  ne  siiinai.  e.  Wi  vh  the  advent  of 

ation  I",  i;  i  me  .'t  aniia.:-  i.is  in  1  ;'tj!:.i  and  tne  need  to  account  tor 

nilii'  liu.  ;l.  1  or  1  111  .1  1 1.  mro  I  y  iTi  arini^ir  tat  U  !' !.  deter  iiii  r :  at  i.  on  .  tne 

h  '  H  nersi;  to  i::l  'tO'  on  no  its  own  set  ot  coordinates  O  t-  tne 

Ijo'e  '-■:::irii  I  ■  i  1 1...  .;i  e  data  supplied  t  y  participating 
■;:'d ter  i  oc  .  'Jn  t 1  t  97  2,  ttn.-*  BlH  used  ornlv'  the  same  data 
■OS  f:.ri'.v,  li  t'o  --a  a  complete  ovs^riap  ol  tunctioris  with  the  IPMS 
I-:-'  !,  s  ted ,  birice  !,  9,1’,  pole  positions  based  on  Doppler 

I  I  cris  i;;f  the  f-ta^s  Na--' i  ijat  i  on  Satellite  Svstem  iNNSS.' 
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sat.ellites  have  been  incorporated  into  the  BIH  global 
•solution.  The  methods  used  by  the  BIH  in  the  treatment  of 
data  are  given  by  Feissel  til. 

Research  at  the  Naval  Weapons  Laboratory  mow  the  Naval 
Surface  Weapons  Center)  by  Anderle  and  Beuglass  C2] 
demonstrated  that  it  was  possible  to  use  Doppler 
observations  of  NNSS  satellites  to  compute  pole  positions. 
Doppler  solutions  of  pole  position  have  been  distributed  bv 
the  Dahlgren  Polar  Monitoring  Service  since  1969.  The  pole 
positions  are  a  byproduct  ot  the  orbit  computation  process. 
Hence,  when  the  responsibility  of  computing  the  NNSS 
satellite  orbits  was  transferred  to  the  Defense  Mapping 
Agency  (DMA)  in  1975,  the  derivation  of  pole  position  was 
continued  by  DMA.  The  DMAHfC  Polar  Monitoring  Service 
(DPMS)  reports  are  distributed  by  the  Hydrographic/ 
Topographic  Center  of  DMA  to  users  on  a  weekly  or  monthly 
b<a3i  s . 

Project  merit  C3J  was  a  program  of  international 
col  lataoration  to  monitor  Earth  rotation  and  mtercompare 
techniques  of  odservatron  and  analysis.  The  MERIT  Main 
Campaign  of  observations  was  held  during  the  period 
September  1,  1933,  to  October  31,  1984  ana  included  a 
variety  of  te'ihniques  for  determini.ng  polar  coordinates. 


2  .  (J  COMPUTAriar  iAL  ME  I  HODS  FOR  OOPfd.ER  FA 

TT>e  pole  determination  metnoci  utilized  bv  DMA  is  the 
method  adopted  bv  the  Naval  Surface  Weapons  Center  iNSWC;  in 
August  1971.  The  brief  aescription  of  the  method  that 
follows  13  taker!  from  the  detailed  descr  i.  pt  i  cn  of  the 
odsei  ■ /at  1  Qfia  i  procedures  and  the  data  reouction  techniques 
gi  ven  i..)  v  Ar.dcrl  j  f.  4  1 , 

Ctjfif:).;.  5-"“  oc  ser ','3  t  . i  iS  a:'€'  mans’  ciai  t  v  r  •/  a  I’lecworK  of 
app-''cv.' 1  ma  tel  .‘itir  I.  c  yi  vie  tr  i  n>.i  stai.ions  controlled  by 

CNri  or',i::l  a  i  i  wtir- k  4  l.i.-.,  ti 'ac  t  i  i  iC  star  i. -ns  ::or'it'"C'i  1  ed  oy 

tne  riyvv  :4Btr  on  ar  1 1  br' oup  ,  i.  lobsei'- vat  i  <jns~  taten  in  a 

r '.)ur  per  j  .'1  are  p.  ocassier.J  wi  th  the  CELEbT  computer 
Dr-ijpr  a.n  to'l.  e  i,  ea  ■=  i;.  — sguares  sc".  ution  i  =■  cPtained  which 
inr 'i.ides  the  si:  i-iO’ is  tan  t  s  of  t;rbiral  i  ntegrat  i  ori .  a  drag 

S'j.al  f -rtc  i:.  or-  tor  eacfi  day,  a  fr'eguency  an>i  troioospheri  c 

r  -vi,'  vtioi'i  scoi  r;:|  f-n.:tor  -Or  eacT:  r'iass,  the  ti’io  cofliponents 
•  !  :-;i  ■  P  'll  i:>;s'  t  :  'Ori  ,  ftnd  t're  'Otirdi  naLes  or  any  mobile 

.4, I  rri  -- 1  till'.  Tre  1  o  .isc —squares  solution  s  based  on 

’i  ■;■!’-•  h’t'-e'-i  tne  cDsr^fr  v  a  f-;  i  tin  3  ana  C'Snoutferd  data, 

•t  ' .  •-•'r:' . t.  ‘  r.rren  I  ctei.;  saieilite  orbit.  Ine 

.1  '  !  ‘1  'I'.'rv ; t ;  rr-.  r  t  le  e  :.iu a. 1 1  ons  -'if  motioi':  come  trom  the 
,  :  a,r  ri,:  '  s  1.  i.. .  The  integration  scheme  is  a 

I  .rtirr  ii  fi',..,- ,  w  1  t  :  a  OHO— m  i  nuts  Htcpsize.  Tne 

:  n.’iii:: r  a  !  ill  s'l,  1  lii:;  i.  .Kies  tlie  taitn's  qr"  at  1 1  at  1  c.-na  i  field, 

<  tii-i  -i.'f  H' :  ■  tir  .Hi,;  S':-. '  jr  ran  ration  pressura,  I  uni -solar 


gravity  perturbations,  and  solid  Earth  tidal  -forces, 
integration  is  done  in  the  true-o-f-date  system. 


The 


The  accuracy  o+'  the  orbit  determination  process  is 
re-flected  b/  tvgo  statistics:  the  discontinuity  in  the 
eph&meris  at  the  ends  o-^  each  2-dav  span  and  the  root  mean 
square  o-f  the  residuals  of  the  least-squares  adjustment 
process.  The  current  accuracy  of  tne  ephemerides  is 
approximately  three  meters  according  to  Murphy  and  Fell  C6]. 
This  value  agrees  with  the  previous  values  given  by  6owman 
and  Leroy  C73  and  Wooden  CS3. 


During  the  MERIT  Campaign,  four  NNSS  satellites  were 
used  at  DliAHTC:  three  “Oscar — type"  satellites  (DMA 
60=1967 92A,  DMA  68=1970-67A,  and  DMA  59=1967-4aA)  and  one 
"IMova-type''  satellite  (DMA  105=19B1-44A>  .  The  Nova 
sivtellite  is  equipped  with  sensors  and  thrusters  which 
compensate  for  non-conservative  forces  in  the  direction  of 
the  velocity  vector.  Thus,  polar  motion  results  from  the 
Nova  satellite  are  largely  unaffected  by  variations  in  solar 
acti  •-'!  tv. 
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AHALYSIS  OF  paPFLER  POLAR  MOTION  VALUES 


The  most  important  components  of  polar  motion  are  terms 
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3.1.  COMPARISON  CJF  OSCAR -Tv PE  AND  NQVrt-TVPE  NHSS  SATELLITES 

The?  di  f  ierences  between  the  bi--dailv  values  tor  each 
individual  NNSS  satellite  and  the  BIH  Circular-D  values  tor 
the  MERIT  period  are  shown  in  Figure  1.  Unfortunately,  only 
the  Nova  satellite  (DMA  lOfi)  has  observations  -for  the  entire 
MERIT  tirnespan.  Although  the  v.  and  v  curves  tor  this 

satellite  have  di+terent  biases,  they  stiow  similar  behavior 
with  a  time  offset  of  approx i matel y  30  days.  Tne  x 

components  of  Oscar  satellites,  DMA  59  and  DMA  bO,  tend  to 
follow  the  Nova  satellite  curve,  while  DMA  *>8  appears  to  have 
periodic  oscillations  about  a  constant  offset.  The  y 

components  of  the  Oscar  satellites  show  less  variation  than 
the  ;!  components.  The  scatter  associated  with  the  Nova 
satellite  values  is  much  smaller  than  the  scatter  in  the 
Oscar  satellite?  values.  The  statistics  associated  with 
these  comparisons  are  given  in  Table  1.  The  means  and 
standard  deviations  in  the  table  are  in  arc  seconds.  The 

number  of  observations  is  given  in  the  column  labeled  n. 
The  second  set  of  means  and  standard  deviations  in  Table  1 
taken  from  the  stud/  by  Colquitt,  Stein,  and  Anderle  C93. 


TABLE  1.  STATISTICS  FOR  THE  OSCAR  -  NOVA  COMPARISONS 
SERIES  I  X-MEAN  SIGMA-X  V-MEAN  SIGMA-Y  N  I  X-MEAN  SIGMA-  X  Y-MEAN  SIGMA-V 
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■•ual  sate!.,  j, 'e-,.  r  .j‘-..'c !.  al  data  sets  oi  5--dav  means  were 
from  th  a  i .  -  ■  ,i  i,  ./r  b  i -da  i  j.  values  for-  the  Oscar 

.  fed  aid  f'./t  i.c;.  r-'i/.a  satellite  to  make  ccmcari  sens  on 
arriar.:!  priLfr-  ir.'  -  .■Dr,  iia/s.  these?  two  data  ser.s  .are 

?tc'i,l  tiO'/A  <r.  .'i  irSCi-.R.  Using  tiiese  si>ec  ai.  oata  sets, 
anci  the  annual,  and  ( rianij.!. ei'  motions  were  removed 
or.)  ser  vat:  j.  I  •;)■•••  Equation  i.  Figure  2 

the  Tj.j.  )•  for  r^nc  :■  3  r -it  we -an  the  residuals  of  the  fits, 
df/si.  r.li  fr.-'.^nti.-v  occ-i.r  between  200  and  davs  into 

■’.iT  i '-iii.'Si.)  ?h .  Ti-. .?  st-atisti-cs  associated  wi  fn  these 
I  i  L  ■'  ,11,  .r'i‘  in  t  1?  next  section. 
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FIGURE  1.  DIFFERENCES  BETWEEN  NNSS  AND  BIN  POLE  VALUES 


X  -  component 


y  -  component 


FIGUEIE  2 
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ia  I  emovf  thK  annual  and  Chandler  information.  The  five 
i  I  oe  par  amet  er  f:.  of  tht'  fite  were  the  meaut  position  of  the 
pole  and  the  amplitudes  of  tlte  sine  and  cosine  terms  of  the 
annual  and  Chiandler  motions.  Table  2  gives  the  results  of 
these  fits  for  each  data  type.  The  entries  of  Table  2  are 
e;;presse?d  in  arc  seconds.  The  BIH  series  is  the  Circular-t> 
final  values.  SLR  is  the  satellite  laser  ranging  series  of 
the  University  of  Texas.  VLBI  is  the  Very  Long  Baseline 
Interferometry  series  of  the  National  Geodetic  Survey.  The 
SLR  and  VLBI  data  were  taken  from  a  magnetic  tape  provided 
by  the  Naval  Observatory.-  DMA  is  the  normal  combined 
solution  of  weighted  means  of  individual  satellites 
published  by  DMA  on  5-day  intervals  in  the  DPMS  reports. 
NOVA  and  OSCAR  are  the  special  5-day  means  discussed 
previ ousl y . 


TABLE  2.  AMF'Ll  TUBES  OF  THE  FITTED  PARAMETERS 
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-iqure  3 

shows 

the 

residuals 

of  the  fits  whose 

f i tted 

parameters 

are 

qi  ven 

i  n 

Table 

^  • 

The 

smooth  line  in 

each 

plot  represents  a  7-point 
being  plotted.  As  is 
geophysical  information  is 
the  major  peaks  agree. 


moving 
evi dent 
simi lar 
All 


residuals  exhibit  a  semiannual 
residuals  show  an  extra  hump  at 
does  not  exhibit  a  semiannual 
signature.  The  differences  in 
and  300  days.  BIH  and  VLBI 


average  of  the  quantity 
from  the  figure,  the 
for  all  data  types,  i.e.  , 
the  plots  of  the  V  value 
period  although  the  OSCAR 
200  days.  The  X  component 
period,  but  has  a  distinct 
the  plots  occur  between  200 
are  relatively  flat  in  that 


interval,  whereas 


BLR  has  its  peak 
days.  Again  the 


all  the  other  data  sets  have  a  small  peak, 
at  275  days  and  NOVA  has  its  peak  at  230 
OSCAR  residuals  are  most  different  with  a 
large  peak  at  215  days.  Paquet,  Djurovic,  and  Techy  tlOT 
also  note  the  similarity  of  the  geophysical  information  in 
their  residual  plots. 


Vo  further  differentiate  subtle  differences  between  the 
Earth  orientation  parameter  series  shown,  in  Figure  3,  the 
residuals  of  the  fits  were  differenced.  The  results  of  this 
process  are  shown  in  Figures  4  through  9.  Note  that  Figure 
2  presents  results  of  this  process  for  the  NOVA  and  OSCAR 
ser 1 es. 


4. 0  PRELIMINARY  EVALUATION  OF  POLE  POSITIONS  COMPUTED  USING 
WORLD  GEODETIC  SYSTEM  1934 

Biases  between  polar  motion  solutions  using  different 
NNSS  satellites  have  been  shown  to  result  from  uncertainties 
xn  knowledge  of  the  Earth's  gravity  field  C113.  The 
question  arises  as  to  whether  the  development  of  the  World 
Geodetic  Svstem  1984  tWGS  84)  and  its  associated  Earth 
Gravi tational  Model  (EGM) ,  as  discussed  by  Decker  C12]  and 
White  r.  13  3,  cari  improve  the  pole  position  solutions.  As 
part  of  ifceststef  ephemeris  generati on } usi nq  WGS  84  (see,  for 
e;:ample,  C14]),  a  number  of  NNSS  ephemerides  based  on  WGS  84 
have  been  computed  at  DMAHTC  using  the  CELEST  program.  Of 
particular  interest  was  an  i nterval,  begi nni ng  in  late  May 
and  ending  in  early  July  1985,  within  which  WGS  84  precise 
ephemerides  for  5  NNSS  satellites  were  produced  on  a  regular 
basis.  At  the  same  time,  analogous  precise  ephemerides  in 
the  NWL  9Z2  s'/stem  using  the  NWL  lOE-1  EGM  were  produced  as 
part  of  DMAHTC  routine  production.  This  situation  permitted 
a  preliminary  comparison  of  the  pole  position  solutions 
computed  in  the  two  systems - 


*+  •  1  FROCEOURES  FOR  GENERA  f  I  MG  WGS  S4  BASED  SOLUTIONS 

WGS  84  versiorts  of  CEL£ST  input  files  were  utilised  to 
.oroduce  pracise  ephemerides  of  5  NNSS  satel  1 1  tes>  f  rom  day 
14'9  to  day  154  and  from  day  li>l  to  day  179,  m  1985.  The 
program  input  included  the  WGS  84  EGM  truncated  to  degree 
(n.'  and  order  (mi  41^  as  well  as  the  WGS  84  ellipsoid 
parameters.  The  program  also  utilized  a  file  of  WGS  84 
tracking  station  coordinates  which  were  obtained  by  applying 
a  geometric  tr  insf  ormat  i.  on  to  the  "production"  NWL  9Z2 
c: or j 1 na tes ,  Each  orpin  fit  for  satellites  DMA  59,  DMh  7Y, 
f.'fiA  93,  and  DMA  113  'ric'/ai  was  based  or.  2  aays  of  tracking 

.■i.ii  each  .  t  to'-  satellite;  DMA  J '>5  was  based  on  1 
■lav'  or  dat-t, 

!/’..i.rinu  the  ’•►est  cerioa.  DiTnri  i  C  aiso  produced  precise 
aphamerides  of  trie  s  .^me  5  iINGS  satellites  for  routine 
production  purpcsos.  These  fits  utilized  the  NWl.  lOE— L  EGM 
(n  =  28,  m-’-27;  .ari’i  the  NWU  97  2  tracking  station  coordinate 
set.  Data  editing  occurred  i ndependent 1 v  in  the  analogous 
CFJ..E3i  solutions  in  the  two  systems  ii.e.,  no  attempt  was 
rti-rie  to  enforce  lomiaon  oaia  between  the  soJutions,'. 

The  piOie  coordinates  output  in  the  two  svstems  were 
crotparnc)  <o  c  I.  A'  iirocuj  ar~D  valuer;,  of  and  y.  Dailv  values 
wer  1.  trr:,ro  tins  ui:'  5 -d  ay  values  by  cuoic  sol  me 

i.  1 1 1  e- pi ,  I  rif  ,1,  on  auf'i  d  !  f  f  er  erices  were  formed  between  the 
i  n  ter  pc  I  ateo  ■-aL  •-vs  and  the  Doppler  ••leri  ved  values.  These 
d f  f .  :r  .aficos  wei  e  cjofjuted  in  the  ser'Se  Doppler"  minus  EUH. 
'h-:-  pole  P'PSition  solutions  which  resulted  from  the 


Tattle  3 :  Statistics  of  the  Dlffeiences  Berveen  the  Doppler-Dexived 
Pole  Cooidinates  end  BIH  Pole  Cktoidinaies 
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Table  4:  Statistics  of  Pole  Cooidtoate  Diffexences  Oxov^  By 

mss  Types 
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OSCARS 
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Notes:  Z  •  mean  difference 
o  «  stendaid  deviation 
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fi 


33b<i>llite  DMA  105  fits  (in  both  systems^  on  day  175  were 
outliers  and  were  rejected  from  further  study. 


PRELIMINARY  RESULTS  OF  WGS  84  -  BASED  POLE  SOLUTIONS 

Table  3  presents  statistics  of  the  differences  between 
the  BIH  and  the  Doppler-derived  pole  coordinates  on  a 
sate! 1 i te-by-satel 1 1 tes  basis.  An  examination  of  this  table 
pave  no  conclusive  evidence  that  use  of  WGS  84 
systematically  improved  the  agreement  between  the 
Doppler-derived  values  and  the  BIH  values.  However,  when 
the  satellites  were  grouped  by  type  (i.e. ,  Nova  and  Oscar), 
as  in  Table  4,  a  possible  correlation  emerged.  While  the 
pole  positions  determined  using  the  Oscar  satellites  and  WGS 
84  showed  no  improvement  in  agreement  with  the  BIH  values, 
pole  positions  determined  using  the  Nova  satellites  and  WGS 
84  did  show  better  agreement  with  BIH.  Additionally,  there 
was  less  scatter  in  the  Nova  differences  when  WGS  84  was 
utilised,  as  evidenced  by  the  smaller  standard  deviations. 


5,0  CONCLUSIONS 

The  polar  motion  derived  from  Doppler  tracking  of  the 
NMSS  satellites  has  provided  and  continues  to  provide 
vaJ  Liable  Earth  or  i  entat  i  on  inf  ormation.  When  the  annual  and 
Chandler  components  of  pal  ir  motion  are  removed  from  the 
obser vat  1  ona I  data  of  e^ach  Earth  orientation  series,  the 
same  geophysical  information  remains.  Nova-type  NNS3 
satellites  provide  more  precise  estimates  of  pole  position 
than  Oscai — type  NNSS  satellites.  Preliminary  analysis  shows 
that  the  use  of  Wi3S  64  may  improve  the  accuracy  of  the  pole 
positicns  derived  from  the  Nova  satellite  orbit  fits.  but 
that  I  s  use  for  Oscar  satellites  does  not  lead  to  improved 
acc'jrac  1  e -1 ,  Oowever  ,  more  data  need  to  be  processed  to 
suh'stan  i  at  tr.ese  pr  c  1. 1.  mi  narv  results. 
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